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Q Abstract 

t£h We consider the billiard dynamics in a non-compact set of M. d that is con- 

structed as a bi-infinite chain of translated copies of the same d-dimensional 
polytope. A random configuration of semi-dispersing scatterers is placed in 
each copy. The ensemble of dynamical systems thus defined, one for each 

£NJ global realization of the scatterers, is called quenched random Lorentz tube. 

Under some fairly general conditions, we prove that every system in the en- 
semble is hyperbolic and almost every system is recurrent, ergodic, and enjoys 
some higher chaotic properties. 

\o 
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o 

1 Introduction 

> 

^ A (i-dimensional Lorentz tube (LT) is a Lorentz gas, in Euclidean <i-space, that is 

confined to a subset T which is infinitely extended in one dimension. 

As a prototype, think of an infinite square-section cylinder in IR 3 , in whose in- 
terior a countable number of convex scatterers are placed approximately with the 
same density (see Fig. [I]). A material point travels inertially in the free region of T, 
until it collides with either a scatterer or the boundary of T (from now on, the latter 
will be referred to as a scatterer as well). Assuming the scatterer to be infinitely 
massive, the collision is totally elastic, i.e., the outgoing velocity v + is derived from 
the incoming velocity v~ by reversing the normal component of v~ relative to the 
plane of collision. 
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In the terminology of dynamical systems, a system like this is an extended semi- 
dispersing billiard. The term 'extended' refers to the fact that the configuration 
space is not compact and the relevant physical measure on it is infinite. Also, it is a 
semi-dispersing billiard because the particle undergoes a billiard-like dynamics with 
bouncing walls that are either flat or convex, as seen from the particle. (We invite 
the reader to avoid confusion with billiards that are semi- dispersing because the 
sectional curvature of the scatterers can be either positive or zero; e.g., 3D billiards 
with cylindrical scatterers.) 

It is a celebrated fact that semi-dispersing billiards give rise to chaotic dynamics 
[CM] — using the term 'chaotic' in a very lax sense here — and so it seems sound to 
use models like this to study the motion of small particles (e.g., electrons) in thin 
wires, in whose interior a configuration of obstacles (e.g., atomic nuclei) makes the 
motion chaotic. (This was more or less Lorentz's original motivation |Loj ; cf. also 
[KFl IAACGI ILWWZl IH&alj and references therein.) 

This note is a follow-up to an article that three of the present authors have 
published recently |CLS] . where these types of systems are studied in two dimensions. 
We refer the reader to the introduction of that paper for a better description of the 
physical and mathematical motivations behind this research. 

Here we just outline our main result and its consequences in terms of the stochas- 
tic properties of the dynamics we consider. Our chief interest, as far as this note 
is concerned, is in the recurrence of these types of systems. Recurrence is the most 
basic property one needs to establish of extended systems in order to study their 
chaotic properties. (It is hard to claim that a certain dynamics "randomizes" the 
state of the system — more precisely, decorrelates it from its initial condition — if a 
non-negliglible part of the phase space is made up of trajectories that escape to 
infinity, thus giving no asymptotic contribution to the state of the system in any 
given compact region.) 

In fact, for our LTs, we will see below that recurrence is a sufficient condition 
for a number of stronger ergodic properties. 

For these effectively one-dimensional systems one would expect recurrence to be 
a typical property. To make this point, we define a fairly large and representative 
measured family (in the language of statistical mechanics, an ensemble) of LTs and 
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ask if the typical element is recurrent, in the sense of Poincare (which coincides with 
the intuitive meaning of the word here). 

The family is defined roughly as follows: The tube T is made up of a countable 
number of congruent d-dimensional polytopes (henceforth cells) having two parallel 
and congruent facets (henceforth gates), whereby each cell is attached to its two 
adjacent cells (Fig. [l] is an example of this). In each cell we put a random configu- 
ration of convex scatterers, according to a very general probability law. Each global 
realization of scatterers defines a different LT. This type of structure, in which one 
randomly chooses a dynamical system and then follows its deterministic dynamics, is 
called 'quenched random dynamical system'. So what we have is a quenched random 
LT. 

If the system verifies some geometric conditions (most of which are rather gen- 
eral, some less), we prove that almost surely, in the sense of the probability, the 
LT is recurrent. Since it can be proved that a recurrent LT is also ergodic, and 
that suitable first-return maps are strongly chaotic (at least K- mixing |AA] ) . an 
important corollary of our work is that the typical LT in our family is chaotic. 

As already mentioned, the main contribution of this paper is an extension of the 
results of [CLSJ to dimension d > 2. It is known by experts in the field that semi- 
dispersing billiards in dimension three and higher present specific subtleties and 
difficulties. Therefore, we made an effort to detail the parts of the proofs that deal 
with such difficulties, while giving a looser exposition of the remaning arguments (as 
they can be found elsewhere as well). The paper is organized as follows: Section |2j 
which should be accessible to the reader with a minimal background in dynamical 
systems, contains the precise formulation of the results. Section [3] gives an outline 
of the main proof, partly referring to previous work by some of the present authors. 
In Section [4], which is the most technical part of the paper, we give precise proofs 
for the arguments that are specific to the systems at hand. 

Acknowledgments. We thank Gianluigi Del Magno and Domokos Szasz for useful 
discussions. This work was partially supported by the FIRB-"Futuro in Ricerca" 
Project RBFR08UH60 (MIUR, Italy). 

2 Mathematical formulation of the results 

Consider a closed <i-polytope C G M d that has two parallel and congruent facets. 
Denoting said facets G 1 and G 2 , call r the translation of M. d that takes G 1 into G 2 , 
and define, for n G Z, C n := r n (Co). Each C n is called a cell and T := \J n& iC n is 
called the tube, see Fig. [2j 

For each n, a family of closed, pairwise disjoint, piecewise smooth, convex sets 
O n)i C C n (i = 1, . . . , N) is given. We refer to this family as the local configuration 
of scatterers in the cell C n (note that some O n ^ might be empty, so different cells 
might have a different number of scatterers). This configuration is random, in the 
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Figure 2: Assembling an LT as a chain of cells with gates. 



sense that each n ^ = n ^(£ n ) is a function of the random parameter £ n G Q, where 
Q is some measure space. The sequence £ := (£ n )nez G which thus describes the 
global configuration of scatterers in the tube T, is a stochastic process obeying the 
probability law EE, whose properties are given later, cf. in particular (Al). 

For each realization £ of the process, we consider the billiard in the table Qt := 
T\ Unez Uili O nt i{£ n ). This is the dynamical system (Qe x S^ -1 , <fr\, me), where 
S^ 1 is the unit sphere in W d and : Qe x S' d_1 — >■ x S^ 1 is the billiard flow, 
whereby (qt,v t ) = <f>\(q,v) represents the position and velocity at time t of a point 
particle with initial conditions {q,v), undergoing free motion in the interior of Qe 
and Fresnel collisions at dQi, i.e., if q t G dQi, then 

v t + = v t - - 2(v t - ■ o qt ) o qt , (2.1) 

where o qt is the inner unit normal to dQi at q t . (Notice that in this Hamiltonian 
system the conservation of energy corresponds to the conservation of speed, which 
is thus conventionally fixed to 1.) Lastly, is the Liouville invariant measure 
which, as is well known, is the product of the Lebesgue measure on Qi and the Haar 
measure on 5 a!_1 . 

We call this system the LT corresponding to the realization £, or simply the 
LT £. As £ ranges in the probability space (f2 z , II), we have a random family, or an 
ensemble, of dynamical systems. This structure is referred to as a 'quenched random 
dynamical system'. As we shall see later, the situation is simplified by the fact that 
these dynamical systems can be reformulated in such a way that they all share the 
same phase space and invariant measure. 

In the remainder, whenever there is no risk of ambiguity, we drop the dependence 
on £ from all the notation. Also, we call universal constant any bound that depends 
on none of the quantities explictly or implicitly involved in the inequality at hand, 
in particular on I. 

We assume the following: 
(Al) IT is ergodic for the left shift a : Q z — > Q z . 

(A2) There exists a universal constant K\ G Z + such that (for all realizations £ G 
fl z ) dO n> i is made up of at most K\ compact, connected, uniformly C 3 (w.r.t. 
n,i) subsets of algebraic varieties (SSAVs), which may intersect only at their 
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borders. These borders, which thus have codimension larger than one, will be 
generically referred to as edges. 

(A3) If q is a smooth point of dQ, let k(g) be the second fundamental form of dQ 
at q. There are two universal constants > k m > such that, for all smooth 
q G dQ, either the SSAV that q belongs to is a piece of a hyperplane or 

< k(q) < k M , 

where the inequalities are meant in the sense of the quadratic forms. 

(A4) There exist universal constants L > 0, K 3 G Z + , and rj G (0,7r/2) such that, 
in each portion of trajectory of length (equivalently, duration) L, there are at 
most K% collisions; and at least one collision with a dispersing (i.e., non-flat) 
part of dQ and such that the angle of incidence (relative to the normal at the 
collision point) is less than 7r/2 — 77. Notice that the above implies the so-called 
finite-horizon condition: the free flight is bounded above. 

(A5) A singular trajectory is a trajectory which has tangential collisions or collisions 
with the edges of dQ (in which case it conventionally ends there). Using this 
terminology, we assume that, for a.e. i and all i,j G {1,2}, there is a non- 
singular trajectory entering Cq through G % and leaving it through GK 

The next and last assumption has to do with the well-known fact that a semi- 
dispersing billiard is a discontinuous (and indeed singular) dynamical system. It will 
be formulated in full mathematical rigor in Section [4j after the necessary definitions 
are given. Here we give a descriptive version which will be quickly understood by 
the "hyperbolic billiardist" . We anticipate, however, that this assumption is verified 
for a reasonable class of perturbations of a periodic LT or in the case in which Q is 
finite, that is, the local configuration in each cell is chosen from a finite number of 
possibilities. 

(A6) There exist a universal constant K4 > such that the Lebesgue measure of 
the ^-neighborhood of each smooth piece of the singularity set does not exceed 
K 4 5. 

We then have: 

Theorem 2.1 Under assumptions (Al)-(A6), the quenched random LT is almost 
surely recurrent, that is, for U-a.e. i G fl z , the LT £ is Poincare recurrent. 

For the sake of completeness, we recall what Poincare recurrence means in our 
context: Given a measurable A C Q x S^ -1 , for m-a.e. (q,v) G A, there is an 
unbounded sequence of times tj such that 0* J G A. 
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Remark 2.2 Although the theorem is valid in any dimension d > 2, the assump- 
tions given earlier were designed for the case d > 3: in the two-dimensional case, as 
presented in |CLS] . the hypotheses are substantially weaker. 



Theorem 2.1 has deep implications. For a fixed £, let D be a finite union of the 
SSAVs that make up dQ; cf. (A2). (So, for example, D could be one smooth piece 
of some dO n ^, or the whole of it.) Denote 



V:= {(q,v) e D x S^ 1 \vo q > 0} , (2.2) 

that is, V is the submanifold in phase space corresponding to the post-collisional 
position-velocity pairs (henceforth line elements) based in D. It is apparent that, if 
the LT t is recurrent, the first-return map onto T> is well-defined almost everywhere 
w.r.t. the natural measure on T> (see below): we call this map T-p. The Liouville 
measure m induces a Tp-invariant measure on P, which we denote fi-p. It is well 
known that dfjtx>(q,v) = (v ■ o q )dqdv, where dq is the volume element in dQ and dv 
is the volume, or Haar, element in 5* d_1 [CMJ. It is a consequence of (A2) that D 
has a finite volume so, upon normalization, we may assume that /xx>(X>) = 1. 

To avoid misunderstandings, let us recall that both dynamical systems (Q x 
S^ -1 ,<p t ,m) and (V,Tx>, fiv) depend on the choice of £ — the subscript has been 
removed only to lighten the notation. 

An important result is the following: 

Theorem 2.3 If the LT i is recurrent, then (Q x S ,d_1 ,0*,TO) is ergodic. Also, 
for any choice of D as described above, (T>,Tx>, Hx>) is K-mixing (thus mixing and 
ergodic). 



Combining Theorems 2.1 and 2.3, we obtain: 



Corollary 2.4 (Q x S d 1 ,<p t ,m) is ergodic and (V, Tt>, Hv) "is K-mixing for U-a.e. 
choice of £ e fF. 

Since generating an example of an LT that verifies (Al)-(A6) may not be imme- 
diate, we present one in Appendix |A| 



3 Flow of the proofs 



The proofs of Theorems |2.1| and |2.3| follows exactly the same strategy as the corre- 
sponding proofs in |CLS|. IL~2j ILlj. For the convenience of the reader, though, we are 



going to outline them in this section, with particular regard to Theorem 2.1 Some 
of the intermediate results will present complications due to the higher-dimensional 
setting. We will explain how to prove these results in closer detail in Section |4j 
The reader is warned that there are a few changes of notation compared to 

piimicLs]. 
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The first step of the proof of Theorem 2.1 consists in showing the hyperbolic 
properties of each dynamical system in the ensemble. So, for the time being, we 
fix I G and describe the LT I by means of a certain Poincare map which we 
introduce momentarily. For n G Z, denote by {-B n ,j}j=i the collection of all the 
dispersing pieces of boundary in the cell C n (it is important that B n j be a whole 
SSAV among those mentioned in (A2), and not just a portion of it). From our 
hypotheses, K n is bounded above by the universal constant KiN. 

Henceforth, we will indicate the index (n,j) with the symbol a, and the space 
of all such indices with A. Evidently, A is countable. For crei, define 

M a := {(q,v) eB a x S^ 1 | v ■ o q > e} , (3.1) 

where e := cos(7r/2 — rj) and 77 is the universal constant that appears in (A4). (Again, 
o q is the inner unit normal to dQ at the point q.) 

All the above definitions clearly depend on I: let us now reinstate this dependance 
in the notation and denote M.i := {Ja^-Ma- By (A4), M.i is a global cross- 
section for the flow <j)\. We call Ti and fi£, respectively, its Poincare map and the 
invariant measure induced by the Liouville measure m f on M< (of course, up to a 
constant factor, fj,£ has the same density as the measure jxx> introduced in Section [2]). 
Each dynamical system thus defined possesses some basic hyperbolic and ergodic 
properties which we now outline, in a rather undetailed way. 

Theorem 3.1 The following holds for the dynamical system (Aie, Te, fie): 

(a) The system is uniformly hyperbolic w.r.t. the natural metric in Aig. 

(b) There is a hyperbolic structure, in the sense that there exists local stable and 
unstable manifolds (LSUMs) at a.e. point of A4(. Also, the two correspond- 
ing (invariant) foliations, when measured with a Lebesgue- equivalent (d — 1)- 
dimensional measure, are absolutely continuous relative to fig. 

(c) Local ergodicity holds. This means that, for each a, a. a. pairs of points in A4 a 
are connected by a chain of alternating LSUMs that intersect transver sally. 
The intersection points can be chosen out of a predetermined full-measure sub- 
set of .Ma- 
in the second step of the proof we represent the LT I in yet another way, which 

will be more convenient later on. What we do is, we introduce a different cross 
section for the same flow. For n G Z and j G {1,2}, denote by G J n := r n (G^) the 
two gates to the cell C n and by Oj be the inner normal to G J n , relative to C n (thus 
02 = — 01). Consider 

^:={(q,v)eGixS d - 1 \v- Oj >0}, (3.2) 

that is, the collection of (almost) all line elements entering C n from the "left" or 
from the "right", depending on j. The global cross section that we use this time is 
N := Unez U,-=i wm le the corresponding Poincare map we denote TVrf- 
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The gates are sometimes referred to as transparent walls, because, in the theory 
of billiards, the corresponding map has virtually the same properties as an ordinary 
billiard map, such as or Tg. In particular, it preserves a measure \xj^ that has 
the same functional form as the measures fiv and fig. 

We end up with the triple (Af, T\frf> Aw)- Notice that neither the phase space 
nor the measure depend on I, which is precisely what makes this dynamical system 
more convenient than the one previously introduced (and gives further justification 
as to why the whole ensemble is called 'quenched random dynamical system': we 
have a family of maps that are defined on the same space and preserve the same 
measure) . 

At this point, one might ask why the map Ti was defined at all. The reason is, 
we needed to prove Theorem 3.1 first, in order to obtain the corresponding results 
for (A/ - , Tfj-.j, fijy). In fact, it is not hard to verify that the latter system inherits the 
hyperbolic structure of the former: One constructs the local stable manifolds (LSMs) 
and local unstable manifolds (LUMs) of (Af, Tj^-.j, u.jy) as push-forwards, respectively 
pull-backs, of the LSMs and LUMs of (A4g,Ti, /i^) |L2j . This is possible because the 
first system has fewer singular trajectories than the second so, for example, when 
a LSM of the second system is pushed forward by the flow, no cuts occur at all 
positive times — ensuring that the defined push- forward is indeed a LSM. 

It is also rather easy to check that uniform hyperbolicity is maintained. The 
result that we are mostly interested in, however, is the analog of Theorem 3.1 (c): 



Theorem 3.2 The dynamical system (Af, Tj^-i, fij^) is locally ergodic in the follow- 
ing sense: For any nfZ and j G {1,2}, a. a. pairs of points in Ml are connected by 
a chain of alternating LSUMs that intersect transver sally. The intersection points 
can be chosen out of a predetermined full-measure subset of M^. 



In the third step of the proof we use the so-called 'point of view of the particle' 
(PVP). It consists of a finite-measure dynamical system that, together with a suit- 
able observable, describes the dynamics of all the orbits in all the realizations of 
the LT. The idea is that, instead of following a given orbit from one cell to another, 
with every iteration of the dynamics we shift the LT in the direction opposite to the 
orbit's displacement, so that the point always lands in the same cell (conventionally 
Co). We briefly outline the construction of this dynamical system, referring the 
reader to |CLS1 IL~2] for more detailed explanations. 

Let ATq := Aq 1 UJVq be the cross-section corresponding to the gates of Co, and 
/i the normalized billiard measure on it. For a given u G Q, determining the local 
configuration in Co, define a map : Ao — > Ao as follows. Trace the forward 
trajectory of (q, v) G Ao until it crosses G 1 or G 2 for the first time (almost all 
trajectories do). This occurs at a point qi with velocity v\. If, for e G {—1, +1}, C e 
is the cell that the particle enters upon leaving Co, define 



Ru(q,v) := (r e (gi),wi). 



(3.3) 
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Clearly R w (q,v) G No, and preserves /i f° r every u. Next define the so-called 
exit function e : No x — > {— 1, +1} via the formula e(q, v; u) := e. From now on, 
we indicate line elements (q,v) with the letter x. 
The PVP system (E, F, A) is defined by: 

• E := N x 

• F(x,£) := (Re (x), a e( ^ x,t °\£)), which defines a map E — E. Here £o is the 
th component of £ and a is the left shift on f2 z , introduced in (Al) (therefore 
^W = KW- with C:=C +e ). 

• A := no x n. Clearly A(E) = 1. Also, using that F is invertible, R^ preserves 
Ho, and a preserves II, it can be seen that F preserves A. 

Now, tolerating the abuse of notation whereby e(x,£) = e(x,£o), let us think 
of the exit function e as an integer- valued observable of the dynamical system just 
defined. We are interested in its cocycle (namely, Birkhoff sum) {SVijneN, given by 
Sq(x, £) = and 

n-l 

S n (x,£) :=^(eoF fe )(x,£). (3.4) 

k=0 

A discrete cocycle, such as {S n }, is said to be recurrent if, for a.e. (x, £), there exists 
a subsequence {rij} such that S nj (x,£) = 0, for all j. For one- dimensional (i.e., 
Z-valued) cocycles, a sufficient condition for recurrence has long been known (see, 

e.g., my 



Proposition 3.3 If (E, F, A) is ergodic, and e : E — > Z is integrable with J s e = 0, 
then the corresponding cocycle is recurrent. 

(A beautiful <i-dimensional version of this result was given by Schmidt [S] — see 
also a generalization in the Appendix of |L3j .) It is not hard to verify that the 



recurrence of {S n } implies Theorem 2.1 In fact, let us call a global configuration 
£ typical if, for all k G Z and /iQ-a.a. x G No, {S n (x, cr fc (£))} n£ N is recurrent. By 
the recurrence of the cocycle, Fubini's Theorem and the denumerability of Z, Il-a.e. 
£ G f2 z is typical. On the other hand, by the above definition, in a typical LT 
£ almost all orbits come back to the cell where they started. Using the Poincare 
Recurrence Theorem on suitable first-return maps, one easily checks that this is 
equivalent to the Poincare recurrence of (N, Tjy : e, fijy) which, clearly, is the same as 
the recurrence of (Qi x S i ~ l ,(j)\,'mt). In turn, since Nte is a global cross-section, 
that is equivalent to the recurrence of (Nte,Tt, ^) and implies the recurrence of 

Since e is bounded and has zero average (this is clear by time-reversal symmetry), 



in order to derive Theorem 2.1 from Proposition 3.3, what remains to be shown is 



Theorem 3.4 (E, F, A) is ergodic. 
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This is proved precisely as Thm. 4.1 of [CLSJ. The idea is to use Theorem 3.2 
to show that each fiber Mq X {£} of £ is fully contained in one ergodic component 
of the system. In other words, the ergodic decomposiontion is coarser than the 
decomposition into fibers. Assumptions (Al) and (A5) then ensure that a IT-typical 
fiber belongs to the ergodic component of any other typical fiber. 



We finish this section by giving a concise outline of the proof of Theorem 2.3 
following [L~1~1IL~2]. 

If (J\f,Tj\f : £, fij\f) is recurrent, then the first-return map to any Ml is well-defined 



almost everywhere. Theorem |3.2| ensures that it is also ergodic, since a. a. points on 
a LSUM belong to the same ergodic component. This implies that no Ml can be 
split into two invariant sets of Tj^-.j. In other words, the ergodic decomposition of 
Ttf.j is coarser than the partition of M into connected components. But (A5) entails 
that Tj^-.i carries a positive measure of points from Ml into M ] n ±\ and M^ 1 (using 
the plus sign for j = 1 and the minus sign for j = 2). This ensures that there is 
only one ergodic component and (M,Tj^ : £, fi^) is ergodic. The ergodicity of both 
(Q x S^ -1 , <p t , m) and (V, Tx>, /i©) follows immediately. 

As for the A-mixing property of the latter dynamical system, once again we con- 
struct its LSUMs as push-forwards or pull-backs of the LSUMs of, say, (M, TV:£, Aw)- 
Here is where the hypothesis that D is made up of whole smooth boundary compo- 
nents comes into play: no further cuts must occur during the push-forward/pull-back 
process. Once a hyperbolic structure has been established for (V, Tx>, Ht>) one uses 
the general result of Pesin's theory [Pj whereby the system decomposes into a count- 
able number of positive- measure ergodic components, over which a power of the map 
is piecewise A-mixing [KSJ. Since in our case it is easy to prove that is ergodic 
for all k G Z + (the above-defined LSUMs are LSUMs for any power of the map as 
well), this immediately implies that the system is A-mixing. 



4 A few detailed arguments 



We have thus seen that, whenever Theorem 3.1 holds, the desired result follows 



by fairly general arguments. This final section — which is the original core of this 



note — is devoted to demonstrating that, under the stated assumptions, Theorem |3.1 
does hold for LTs in dimension 3 and higher. 

The problem with hyperbolic billiards in d > 3, as is common knowledge in 
the field and perhaps not so common elsewhere, is that the so-called fundamental 



theorem (namely, local ergodicity; cf. Theorem 3.1| fcj) is not known for general 



semi-dispersing tables, even in finite measure. The misinformation is due to the fact 
that incorrect proofs of said theorem were believed valid until recently, when Balint, 
Chernov, Szasz and Toth |BCST] pointed out the mistake. In the same paper, these 
authors recover the proof for the case of algebraic Sinai billiards, i.e., dispersing 
billiards on the torus with a finite number of scatterers, whose boundaries are made 
up of a finite number or compact pieces of algebraic varieties. To our knowledge — if 
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we exclude generic results that so far can claim no specific examples |BBTj — these 
are essentially the only semi-dispersing tables, in d > 3, for which ergodicity is 
known. 

Because the situation for Sinai billiards is less than optimal, our ability to prove 
local ergodicity for our models is also less than optimal. In truth, we simply adapt 
the results of |BCSTj to the framework at hand, much as a previous paper by one 
of us |L1] adapted the classical results on two-dimensional semi- dispersing billiards 
to planar Lorentz gases. 

In what follows, we drop the subscript i from all the notation, denoting our 
dynamical system as (Ai,T,n). To start with, one needs to establish uniform hy- 
perbolicity, namely, statement (a) of Theorem 3.1, because, as we shall see below, 
that is used in the proof of (c). 

In 2D, uniform hyperbolicity descends rather easily from the fact that, after a 
certain time past a dispersing collision (with curvature bounded below), any tra- 
jectory has acquired a sufficient amount of hyperbolicity; namely, if one constructs 
an infinitesimal dispersing beam around the trajectory, the beam has increased its 
dispersion by a large enough factor, cf. |CLSj . 

For d > 3, a complication arises that is sometimes referred to as the problem 
of astigmatism [BJ. It turns out that dispersing beams of trajectories that collide 
with a scatterer almost tangentially (these are usually called grazing beams) do not 
acquire much additional dispersion. This problem is circumvented if one prescribes 
that, a positive percentage of the time, any given trajectory undergoes a dispersing 
collision that has a non-negligible head-on component |BDj . 

Assumption (A4) guarantees this. In fact, Ai is defined as the cross-section of 
all line elements on a dispersing boundary whose outgoing velocity has a sufficiently 
large component along the normal vector to the boundary. Ai is thus a global cross 
section, with corresponding map T. The same assumption also guarantees that, 
after K 3 returns to Ai, a trajectory has traveled at least a distance L, and thus (by 
(A3) as well) has acquired enough hyperbolicity, relative to the so-called orthogonal 
Jacobi metric in tangent space [W] IBP] . But, on Ai, because of the inequality in 
(3.1), this metric is equivalent to the natural Riemannian metric. Therefore T Ks is 
uniformly hyperbolic. But this implies the same for T, since the orthogonal Jacobi 
metric is non-decreasing for tangent vectors corresponding to dispersing beams. 

In order to give a rigorous formulation of assumption (A6) and explain how it 
is used in the proof of Theorem 3.1| f6 ), we need to lay out some facts and a bit of 
extra notation. 

It is common knowledge that billiard maps such as T are discontinuous. If x G Ai 
is the initial condition of a singular trajectory that has a tangential collision or hits 
an edge before the next return to Ai, then quite generally x is a discontinuity point 
of T. We call such x a singular point for the map T. (If x is singular because of a 
tangential collision, it can be seen that the differential of T blows up at x, whence 
the term 'singular'.) 

In our case, given the peculiar choice of Ai, cf. (3.1 ), we must consider as singular 
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also those trajectories which, after a collision, have a velocity v such that 

v ■ o q = e. (4.1) 

The corresponding line element (q, v) is clearly a discontinuity point for the Poincare 
map of M.. 

Let S denote the set of all singular points of T and define S a := S D -M Q . It is a 
well-known and easily derivable fact that S a is decomposed into smooth portions of 
codimension-one manifolds, each of which corresponds to a source of singularity (a 



tangential scattering, an edge, or condition (4.1)) encountered before or at the next 
return to M.\ and to the itinerary of scatterers visited before that. By (A2) and 
(A4), the number of scatterers (and thus number of edges) that can be visited before 
the next return to M. is bounded by a universal constant. Therefore the number of 
smooth pieces that S a comprises is also bounded by a universal constant. 

Moreover, since the LT is algebraic in the sense of (A2), an easy adaptation of the 
results of |BCSTj implies that S a is actually a finite union of SSAVs, whose number is 
universally bounded. (The proof of the algebraicity of the singularity set, in [BCST, 
§5.1], does not use in an essential way that the scatterer configuration is periodic 



there.) Notice that the extra singularities due to (4.1) also give rise to SSAVs 



(Substitute eqn. (5.4) of [BCSTJ with the polynomial equation corresponding to 



(4.1). 



For 5 > 0, define 

S [ * ] := {x e M a | dist(x, S a ) < 5} . (4.2) 

The measures of these neighborhoods play a pivotal role in the proof of the hyper- 
bolic properties of billiards. The considerations in the previous paragraph and the 
results of |BCSTl §5.2] imply that, as 5 — > 0, 

Leb(5W) = 0(S), (4.3) 

where Leb is the Lebesgue measure on Ai (more precisely, the Riemannian volume 
on M. corresponding to the distance dist; notice that /i is absolutely continuous 



w.r.t. Leb.) The implicit constant in the r.h.s. of (4.3) depends in general on £ and 



a: we require the bound to be uniform. More precisely we reformulate: 

(A6) There exists a universal constant > such that, for all sufficiently small 5, 
Leb{S l a ] ) < K A 5. 



By (4.3) it is not hard to generate examples of LTs satisfying (A6). For example, 
one can start with a periodic algebraic LT and then perform a (quenched random) 
algebraic perturbation. By this we mean that the equations of the perturbed scatter- 
ers are polynomials whose coefficients are very close to the corresponding coefficients 
for the unperturbed cell. Another easy example is the case where Q is finite. In 
that case, say that B a is contained in the cell C n . Then, by (A4), S a is completely 
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determined by the local configurations in the cells C^, with n — K 5 < k < n + K 5 
(K 5 being a universal constant). But there are only finitely many possibilities for 



these local configurations, therefore (A6) is implied by (4.3) 



The argument that proves Theorem |3.1(?>j is virtually the same as in Lem. 3.2 



of |L1] . It is based on the old principle that, in good hyperbolic billiards, the 
exponential expansion of the LUMs, or candidates therefor, is the dominant effect, 
compared to the cutting operated by the singularities. Thus, locally along a given 
orbit, one has all the ingredients of Pesin's theory to prove the existence of a LSM. 
The same ingredients then guarantee the absolute continuity of the correspondong 
foliation. 

More specifically, we want to show for /i-a.a. x G M, a constant C = C (x) can 
be found such that 

dist (T~ k x, S U dM) > C k~ 3 , (4.4) 

for all positive integers k. Let us fix x G M a . By the finite-horizon condition, T~ k x 
can only belong to a limited portion of the phase space, namely U/3ev fe M-p, where 
Vfe = Vk{ot) is the index set of all the boundaries that can be visited within time k 
by a trajectory starting in B a . Using (A2) as well, 

#V k <K 6 k, (4.5) 

for some Kq > 0. So we have that the statement 

dist(T- fc x, S U dM) < k~ 3 (4.6) 

is equivalent to the statement 

xeT k ( |J (SpUdM^'A. (4.7) 

\/3ev fc / 



By the invariance of /i, (4.5) and (A6), the measure of the r.h.s. of (4.7) is bounded 
by a constant times k~ 2 , which is a summable series in k. By Borel-Cantelli applied 



to the finite-measure space (M. a , /i), the event (4.7), equivalently (4.6), may happen 
infinitely often in k only for a negligible set of x, whence (4.4). The same reasoning 
applies of course to every a. 

Finally, for the statement (c), we prove our version of the fundamental theorem 
(cf. |L1} Sec. 4]) using the technique of regular coverings, as in [KSSJ or |LW] . This 
technique requires a global argument (i.e., an estimate on objects outside of the 
neighborhood under consideration) in one part only, the so-called tail bound. The 
rest of the proof is local, thus unable to distinguish between a finite- and an infinite- 
measure billiard: all the standard arguments — including the exacting ones where 
one uses that the singularity set and its images via the map are made up of SSAVs 
|BCSTj — apply there. 
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The tail bound is the following statement: For all i 6 A1, there exists a neigh- 
borhood U of x such that 



/u < x e U 



distpi/ s x, 



|J T~ k s\ 

k>m / 



6o(l) 



(4.8) 



as m — > oo. Here distvK s (^, •) is the Riemannian distance along W s (x). (Compare 
(4.8) with the statement of Lem. 4.4 of |L1] . noticing that here we use T~ l and S, 
instead of T and S~ , the latter denoting the singularity set of T~ x .) Once again, 
there is no loss of generality in choosing U C M. a . By the earlier reasoning, the only 
singularities whose images via T~ k can get close to x e Ai a are those in U^ev fe 



Therefore (4.8) descends from the estimate: 



/i < x e M c 



distyt/ s I X 



< 



|J < x e M a 
|J IxeM 



dist^ I x 



distws I T k x 



k>m (3eV k / 

T- k \JsA<t 

/?ev fc / 
U SA <5c\ 



< E Jr' Us 

fc=m+l \ \/3eV fc 
oo 

< const 5 ^ A;A fc . 

fc=m+l 



[<5cA ft 



(4.9) 



In the first inequality we have used the uniform hyperbolicity of T (A < 1 is the 
contraction rate and c is a suitable constant). The third inequality follows from the 
invariance of \i and (A6). 



A Appendix: An example of a 3D Lorentz tube 

In this appendix we present an example of a three-dimensional LT that verifies all 
the assumptions of Section [2] 

We begin by constructing a template cell with its set of scatterers. Then we check 
that the geometric assumptions (A2)-(A4) hold for this cell. Finally, we specify how 
to use that to construct a quenched random LT that verifies all assumptions. 

The cell is a rectangular parallelepiped with a square base of side length 1 and 
with height h > 1. For the sake of the description, we embed this solid in the 
(x, y, ^)-space: the square base lies in the (y, z)-plane, so that the edges of length 
h are in the x-direction. We call the x-direction 'longitudinal' and any direction 
orthogonal to it 'transversal'. 



Lorentz tubes 



15 



Puncture each square facet with a square hole of side g, centered in the center 
of the facet. These holes take the role of the gates G l . Populate the cell with the 
following scatterers: For each of the four edges of length h, consider the cylinder 
of radius p which has that edge as its axis. Choose 1/2 < p < (1— g)/\/2, so 
that the cylinders do not obscure the gates and intersect each other with a positive 
angle (Fig. [3J. The portions of the cylinders that are inside the cell are morally our 
scatterers. However, in order to satisfy assumption (A3), we modify these quarter- 
cylinders by adding a small positive curvature along the longitudinal direction. We 
do so in a way that in every transversal section we still have four quarter-circles 
whose radii satisfy the inequalities given earlier — in other words, we have a so-called 
diamond billiard; cf. Fig. |3J For lack of a better name, we call the resulting solids 
'cigars'. 




Moreover, we insert a 'bulkhead' in the middle of the cell, that is, we add an 
infinitesimally thin scatterer given by the intersection of the plane passing through 
the center of the cell and orthogonal to the vector (1, 1, 1), and the cell itself. We 
punch an off-center, small, polygonal hole through the bulkhead, so that it is possible 
to go from one "chamber" to the other, but no free (i.e., collisionless) trajectories 
exist between the two gates of the cell. Notice that by choosing the parameter h 
large enough, we can ensure that the maximum number of consecutive collisions 
between the flat facets and the bulkhead is 3. 

Next, we argue that this setting verifies (A4). We say that a trajectory that hits 
a scatterer almost tangentially grazes it. So consider a trajectory that grazes a cigar. 
If this trajectory has a sizable transversal component (compared to the longitudinal 
one) then its projection on a transversal section is not far from a grazing orbit in a 
diamond billiard, which implies that the next bounce is not grazing. Let us hence 
consider grazing trajectories that are almost longitudinal. In this case, within a 
certain (bounded above) time, the point will hit the bulkhead in such a way that 
the next collision is necessarily non- grazing and against a cigar, provided h is not 
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too small. On the other hand, since the "transversal angle" at each intersection 
point of two cigars is bounded below, no more that a bounded number of collisions 
can be performed in said time. These arguments will be made more quantitative 
below. 

Finally, the full LT is made up of random algebraic perturbations (in the sense 
of Section [4]) of this local configuration. If the perturbations are i.i.d. in each cell, 
(Al) is obviously verified. If they are sufficiently small, (A6) holds. As for (A5), it 
is not hard to see (by ergodicity of the inner dynamics in a cell, if one will) that this 
assumption is verified as well. 




Figure 4: The longitudinal projection of a trajectory in the (purely cylindri- 
cal) cell, describing a planar trajectory in a diamond billiard. The indicated 
part of the boundary represents a 'zone' (see text). 



At the request of an anonymous referee, we give estimates for the universal 
constants L and K 3 for this LT. We begin by considering a periodic tube made up of 
template cells that — contrary to what we have imposed earlier — have no longitudinal 
curvature. So the dispersing scatterers are truly cylindrical and, if we neglect the 
bulkheads, describe the same diamond in each transversal section (Fig. [4]). Call 7 
the angle at each vertex of the diamond (so 7 = arccos(l/2p)). Fix h > 10, say, and 
g<10-\ 

Denote by v x and v yz , respectively, the longitudinal and trasversal components 
of the velocity v of the material point (thus, \v x \ 2 + \v yz \ 2 = 1). The longitudinal 
projection of a piece of trajectory that does not hit a bulkhead is a (planar) trajectory 
of the diamond billiard. Consider a collision against a cylinder (for the 3D billiard) 
with outgoing velocity v. Denoting by 9 G [0, tt/2) the angle of incidence of v w.r.t. 
the normal at the collision point and by ip the corresponing angle for the projected 
trajectory, we have cos^ = \v yz \ cosip. 
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Let us make a couple of observations on the dynamics in the diamond billiard. 
To start with, we subdivide its boundary in 4 'zones', each zone being defined as 
the set of all the boundary points that are closer to a given vertex than to any other 
(see Fig. [4]). Let M := [77/7], i.e., M is the minimum integer > 77/7. It can be seen 
that any trajectory in the diamond can have at most M consecutive collisions in the 
same zone, and at least one of them will have an angle of incidence (p < (tt — 7)/2 
(to see this, just "linearize" the corners). Also, the time between the first collision 
in a zone and the first collision in another is bounded below by L\ := l/v2 — p (this 
is the semilength of the biggest square inscribed in the diamond). 

Coming to three-dimensional trajectories, we distinguish two cases: the trajec- 
tories with a good transversal component, defined by \v xy \ > 1CT 2 , and the other 
ones, which we call almost longitudinal. Starting with either case, we want to find 
an upper bound for the time before the next "head-on" collision — the parameter 
77 that defines "head-on", cf. (A4), will de determined later. By time-reversibility, 
twice this upper bound will be a good estimate for L. 

If a trajectory keeps a good transversal component during its entire visit in 
a single zone, we know from above that there will be a collision with angle of 
incidence 9 < arccos(10 -2 sin(7/2)) =: tt/2 — rj, in a time less than L 2 : = 2hM 
(2h is an estimate for the horizon of our billiard). If not, we eventually fall in 
the next case, which we consider right away. An almost longitudinal trajectory 
will remain such until it hits a bulkhead, and this occurs necessarily within a time 
L 3 := 3/i(l — lCT 4 ) -1 / 2 . The next collision after that, given the inclination of the 
bulkhead and that h > 10, is against a cylinder, with 9 < n/2 — 77. 

So, L can chosen to be 2(L 2 + L 3 ). During this time, the point can visit at 
most \L/L{\ zones, therefore an upper bound for the number of collisions it can 
have is K 3 := \L/L{\ M + \3L/h] (the first term estimates the collisions against the 
cylinders and the second term the collisions against the flat boundaries). 

We treat the original LT as a perturbation of the one just considered. If the 
radius of curvature of the cigars in the longitudinal direction is, say, larger than 
100/i, multiplying the above estimates by 100 will certainly work. 
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